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Abstract: We introduce an approacli to calculate tine tfiermody- 
namic oxidation and reduction potentials of semiconductors in aque- 
ous solution. By combining a newly-developed ab initio calculation 
for compound formation energy and band alignment with electro- 
chemistry experimental data, this approach can be used to pre- 
dict the stability of almost any compound semiconductor in aque- 
ous solution. 30 photocatalytic semiconductors have been stud- 
ied, and a graph (a simplified Pourbaix diagram) showing their va- 
lence/conduction band levels and oxidation/reduction potentials is 
produced. Based on this graph, we have studied the stabilities and 
trends against the oxidative and reductive photocorrosion for com- 
pound semiconductors. We found that, only metal oxides can be 
thermodynamically stable when used as the n-type photoanodes. 
All the non-oxides are unstable due to easy oxidation by the photo- 
generated holes, but they can be resistant to the reduction by elec- 
trons, thus stable as the p-type photocathodes. 



One key issue in the research of photocatalytic water splitting 
is to search for semiconductor photoelectrodes that can absorb the 
visible light and drive the hydrogen (oxygen) evolution reaction 
using the photo-generated electrons (holes). '^ This requires the 
semiconductors to have proper band alignment relative to the wa- 
ter redox potentials, e.g. the conduction band minimum (CBM) 
of the p-type photocathode should be higher (more negative in po- 
tential) than the water reduction potential H+/H2, and the valence 
band maximum (VBM) of the n-type photoanode lower (more pos- 
itive in potential) than the water oxidation potential O2/H2O, as 
shown in the band alignment plot ([T]l for the Z-scheme water split- 
ting system. ^ Given the importance of the band alignment, it is 
usually taken as a screening condition in the search and design of 
new photocatalytic semiconductors. However, this condition is not 
sufficient, and one important issue, which has attracted increasing 
attention but has not been well studied, is how to evaluate and en- 
hance the stability of semiconductors in the aqueous solution, ^-'^i^' 

Resistance to the photo-induced corrosion (degradation or de- 
composition) under illumination is a critical condition for the pho- 
tocathode and photoanode materials. It is due to this condition 
which makes the photocatalytic water splitting a much more chal- 
lenging problem than the photovoltaic. A compound semiconduc- 
tor MX (e.g. M=Zn, Ga, Ti, X=S, N, O2 in ZnS, GaN, Ti02 respec- 
tively) used as the n-type photoanode may have a VBM lower than 
the O2/H2O oxidation potential, but the photo-generated holes (/;+) 
may oxidize the semiconductor first, rather than the water, making 



the compound MX decomposed through this reaction. 
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MX + zh'^ + solv ^ M-'+ •solv + X {oxidization) 
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Figure 1. A schematic plot of the band alignment of the p-type photo- 
cathode and n-type photoanode semiconductors relative to the water redox 
potentials in the Z-scheme. 0°' shows the oxidation potential of the pho- 
toanode in aqueous solution, and (p" shows the reduction potential of the 
photocathode. 



Similarly the photo-generated electrons (e^) in the p-type photo- 
cathode MX may reduce itself rather than the water through the 
following reaction. 



MX 



-solv^M + X~ •solv {reduction) 



(2) 
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The above reactions define two potentials for the hole and electron 
respectively. When the electron (hole) Fermi energy equals that po- 
tential, the reactions are in equilibrium, i.e. the Gibbs free energy 
change equals zero, and when the photo generated electron (hole) 
quasi Fermi energy is higher (lower) than that potential, the reac- 
tion will occur, hence the semiconductor will be corroded. These 
potentials are called the thermodynamic reduction potential 0'''' and 
oxidization potential (^'" of the semiconductor. 

Since the photoanode (photocathode) is in contact with the aque- 
ous solution, the photo-generated holes in the valence band (elec- 
trons in the conduction band) can oxide (reduce) either the water 
or itself, as shown in |2] Whether the semiconductor is resistant 
to the photocorrosion depends on the alignment of ^"^ relative to 
(|)(02/H20) for the photoanode, and (^'''' relative to ^(H+/H2) for 
the photocathode, as shown in|2] Generally speaking, a semicon- 
ductor is stable with respect to the hole oxidation if its ^"^ is lower 
than either ^(OifRiO) or its VBM, and is stable with respect to the 
electron reduction if its ^'''^ is higher than either (|)(H+/H2) or its 
CBM. 

Experimentally, (^'" and ^''^ can be derived from the Pourbaix 
diagram, I ' but the diagram is not available for many novel 
photoelectrode semiconductors, and it is not easy to measure it. 
Thirty years ago, 0°'^ and (|)'''° of several binary semiconductors 
had been calculated by Gerische r 1 (where the labels pEdecomp 
and n£decomp Correspond to (j)"'^ and 0'''' respectively) and Bard and 
Wrighton. Park and Barber have also calculated the full Pourbaix 
diagram for a few simple semiconductor compounds.— These cal- 
culated <j)°-'^ and (j)'^" are widely cited in literatures to explain the 



corrosion of semiconductors in aqueous solution. 



posed by Lewis and coauthors in their review paper,— the 0°* and 



11,17,18 

1 



As pro- 



(a) photoanode 



(b) photocathode 



/\ 



/\ 



oxidized 



MOj/H.O) 



stable 



oooooooo) VBM 

stable 



r 



stable 



CBM 



stable 



-^(HVHj: 



reduced 



Figure 2. The stability change of the photoanode (a) as its oxidation po- 
tential (|)'" shifts up from below the VBM to above 0(O2/H2O), and of the 
photocathode (b) as its reduction potential ^''' shifts down from above the 
CBM to below 0(H+/H2). 



0''" should be taken into account when choosing the candidate pho- 
toelectrodes, in additional to the optimal band gap and band edge 
alignment. Because the calculation of <!)"■'■ and (|)'''° requires the 
Gibbs free energy of all the reactants and products in the related 
reactions, which are not available for some semiconductors (espe- 
cially for multinary ones), 0"-* and 0''' of only a few simple semi- 
conductors had been reported, and no general trend had been stud- 
ied systematically. Recently new ab initio methods based on the 
density functional theory have been developed to calculate the for- 
mation energies of semiconductors with high accuracy, and also 
more experimental electrochemistry data can been found in the lat- 
est handbooks, ^''•^' thus a more extended study of ^"-'' and ^'''' for 
almost all semiconductors becomes possible. To facilitate the fu- 
ture choice of candidate photocatalytic materials, we here report 
the (^i"* and ^"° for more than 30 semiconductors as well as their 
band alignment relative to the normal hydrogen electrode (NHE) 
potential. Based on these results, we will discuss the trend in the 
oxidative/reductive stabilities of a series of metal oxides, II- VI and 
III-V related semiconductors. 

To start with, we will take ZnS as an example and introduce how 
we calculate its (p"^ and (p''" . The specific exemplification of Eqs.(l- 
2) for ZnS are, 

ZnS + 2h+ +H20^ZnO + S + 2H+ {oxidization) (3) 

ZnS + 2e' +2H+ ^Zn + H2S (reduction) (4) 

which define the (j)"'' and (j)'''' of ZnS respectively, i.e. when the 
chemical potentials (Fermi energies) of the holes and electrons are 
equal to (p"" and (f)''" respectively, the Gibbs free energy changes of 
Eq.(3) and (4) are zero, thus, 

(j)"" = [G{ZnO) + G{S) + 2G(H+) - G{ZnS) - G(H20)]/2eF (5) 

(p'" = -[G{Zn) + G(H2S) - GiZnS) - 2G[H+)\/2eF (6) 

where G(X) stands for the Gibbs free energy of X at the stan- 
dard state, e is the elementary charge and F is the Faraday con- 
stant. To reference the potentials to the hydrogen reduction poten- 
tial (^(H+/H2), the following half reactions can be used: 

H2+2h+ ^2H+ 

2H+ 



-2e- 



(7) 
^ H2 (8) 

Here the chemical potentials of both the electron and hole are equal 
to 0(H+/H2), thus /?+ is equivalent to — e^, and the above two equa- 
tions are equivalent. (^(H+/H2) can be written as 

(p{H+/H2) = [2G{H+)-G{H2)]/2eF (9) 

Now, we can add the reverse of Eq.(7) [exchange the left and right 
hand side] to Eq.(3), then get the sum reaction. 



ZnS + H2O ^ ZnO + S + H2 (oxidization) 



(10) 



and add the reverse of Eq.(8) to Eq.(4), then get 

ZnS + H2 ^ Zn + H2S (reduction) (11) 

Note, to make /i+ in Eq. (3) and (7) and e^ in (4) and (8) equiv- 
alent, it is assumed that their chemical potentials are the same as 
0(H+/H2), so now the Gibbs free energy changes (AG) of Eq. (3) 
and (4) are not zero, and they are equal to AG of the reactions (10) 
and (11) respectively, thus the sign — >■ is used rather than ^. AG of 
the reactions (10) and (11) can be calculated directly, 

AG(IO) = G(ZnO) + G(S) + G(H2) - G{ZnS) - G(H20) (12) 

AG(1 1) = G(Zn) + G{H2S) - G(ZnS) - G(H2) (13) 

From Eq. (5), (6), (9), (1 1) and (12), we have 

(/>"■' = AG(10)/2eF + (^(//+//f2) (14) 

f-' = -KG(n)/2eF + (p(H+ /H2) (15) 

Since the NHE potential is (/i(H+/H2) at pH=0, ip"^ and ip'" rela- 
tive to NHE can be calculated using Eq. (14) and (15). The key 
of the procedure is to find the plausible oxidation and reduction re- 
actions, such as (10) and (11), and get their AG, then the relevant 
^"^ and ^'''' can be calculated. Several plausible reactions can be 
tried, then the lowest (p''" and highest (^"^ should be used as the true 
reduction and oxidization potentials. Thus, in a sense, what we get 
is a simplified Pourbaix diagram, only the potentials limiting the 
photo-corrosions are reported here. In the supplement materials, 
we listed the reactions we considered in determining (^"■* and (^''^ 
for different semiconductors. It should be mentioned that there is no 
guarantee that we have enumerated all the possible reactions, thus 
strictly speaking, the potentials reported here should be considered 
as a higher limit for (|)'''° and a lower limit for 0°'^. Nevertheless, 
in practice, we believe the results we obtained are very close to the 
true oxidation/reduction potentials. 

[3]shows our calculated ^°-'' and ^'''^ relative to NHE, the water re- 
dox potentials (^(02/H20) and (^(H+/H2), and the valence and con- 
duction band edges for five selected classes of semiconductors at 
pH=0. It is well-known that (^(02/H20) and (^(H+/H2) depend on 
the pH value according to Nemstian relation, ^^ i.e. shifting up by - 
0.059 V as the pH increases by 1. The dependence of 0'" and ^'''' on 
the pH value is determined by the specific reactions: when the sum 
reactions such as Eq. (10) and (11) are irrelevant to H+ or OH^, 
their AG is fixed, then ^"-^ and ^''^ shift together with (/)(H+/H2) 
as described by Eq. (14) and (15) and also follow the Nemstian 
relation; when the reactions are relevant to H+ or OH^, the de- 
pendence is different, at variance to the Nemstian relation. Most 
of ^"^ and (^'''' plotted in [3] follow the Nemstian relation, similar 
as (^(H+/H2), and the exceptions include: (i) all oxides have fixed 
(I)"* with respect to pH change, (ii) ^"■'' of CuGaS2, Cu2ZnGeS4 
and Cu2ZnSnS4 shift up by -0.044 V (instead of -0.059 V) as pH 
increases by 1. With these relations clear, the 0'" and ^'''^ relative 
to NHE or (^(02/H20) and (^(H+/H2) at different pH can be plotted 
(as a simplified Pourbaix diagram). 

The valence and conduction band edges in [3] are collected from 
Ref. [Il l for most metal oxides and from Ref. [23,24 1 for non- 
oxides. For the metal oxides, the band edge position relative to 
NHE is reported to depend on the pH value of the solution through 
the same Nemstian relation as for the water redox potentials. For 
the non-oxides, the relation between the band edge position and 
the pH value is so far not clear. ^^ Ref. [ 23 ] derived the band edge 
positions of the group IV, III-V and II- VI semiconductors relative to 
NHE based on the calculated band offsets among semiconductors, 
assuming that the conduction band edge of Si corresponds to its 
electron affinity (about -3.7 eV relative to vacuum level) and NHE 
is -4.44 eV relative to vacuum level. This implicitly assumes that 
the band edge positions for these materials do not depend on the 
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Figure 3. The calculated oxidation potential 0°*^ (red bars) and reduction potential (p'''^ (black bars) relative to the NHE and vacuum level for a series of 
semiconductors in solution at pH=0. The water redox potentials (^(02/H20) and (^(H+/H2) (dashed lines), and the valence (green columns) and conduction 
(blue columns) band edge positions at pH=0 are also plotted. The alignment for these potentials at different pH values (a simplified Pourbaix diagram) can be 
derived according to the relations described in the text, such as the Nemstian relation. 



pH value. The band offsets of almost all group IV, III-V and II- 
VI related semiconductors had been calculated using the ab initio 
methods with reasonable consistence compared to experiments,^ 
and the similar calculation procedures can be used to estimate the 
band edge positions of other semiconductor whose experimental 
values are unavailable. 

Compared to the calculation performed by Gerischer et gZ., i^^ 
Bard and Wrighton, '^ and Park and Barber''* 30 years ago, our 
current work benefits from the progresses in both experiment and 
ab initio calculations. It is the combination of these two progresses 
which allow us to calculate (j)"'' and (p'''^ for almost any given semi- 
conductors, no matter binary or multinary compounds, as shown in 
[3] (i) Experimentally, the Gibbs free energy data in handbooks like 
Ref. [20,21] covers much more chemical compounds with higher 
accuracy compared to those published 30 years ago. (ii) If the ex- 
perimental Gibbs free energy (formation enthalpy plus entropy con- 
tribution) for a compound in the proposed reaction is not reported 
in those handbooks, mostly for some novel compound semiconduc- 
tors, its formation energy can be calculated using modem ab initio 
methods. The Gibbs free energy of that compound is then calcu- 
lated from the formation energy plus the Gibbs free energies of the 
elementary compounds which are all known experimentally. For 
example, in this work, we have calculated the formation energies 
of Ta2N5, TaON, ZnGeP2, Cu2ZnGeS4, Cu2ZnSnS4, etc. using a 
newly developed approach '^ based on the density functional the- 
ory. This approach reduces the formation energy error from 0.25 
eV/atom to less than 0.05 eV/atom for the compound semiconduc- 
tors. '"^ As a result, if a compound in the reactions has 4 atoms per 
formula unit, the errors of the calculated 0°* and tj)'^" are reduced 
from 1 V to 0.2 V (assuming one electron is transfered in the reac- 
tion). Note that the entropy contribution in the formation energy is 
neglected for those semiconductors with no experimental data, but 



for crystals this contribution is small and we estimate it causes er- 
rors less than 0.2 V in (j)"-'' and (j)'''. In the supplement materials, we 
listed how we get the Gibbs free energy of all the solid compounds 
which are not available in Ref. 120,21 1. 

In the following we will discuss the trends in[3]and their influence 
on the design of photocatalytic system. 

(i) For the oxidation potential, most metal oxides have 0°'^ lower 
than (/i(02/Il20) at pH=0, indicating they are resistant to the hole 
oxidation and stable in the solution, but four exceptions are also 
found, CU2O, ZnO, PbO and FeTi03. The reasons for the easier ox- 
idation of CU2O, PbO and FeTiOs are obvious, since the cations in 
these compounds are not at their highest valences and can be further 
oxidized to higher valences (Cu+ to Cu^+, Pb^+ to Pb'*^ and Fe^+ 
to Fe3+).27 Comparing CU2O with Cu+ and CUWO4 with Cu2+, 
we can find that CUWO4 has a much lower 0™ than CU2O, since 
both cations in CUWO4 are in the highest valence state. According 
to this simple trend, all the metal oxide semiconductors at a not- 
highest valence state of cations are likely susceptible to oxidation 
in water, and tend to be unstable (Note C03O4 is an exception, Co 
has higher valence in C02O3, but C02O3 is not as stable as C03O4, 
thus its <j>°-^ is quite low). The easy oxidation of ZnO is also unique 
considering that Zn is at its highest valence state and its band gap is 
quite large with a very low valence band. Thus the stability of ZnO 
cannot be predicted according to its valence state of elements, or its 
valence band position. This is related more to its formation energy 
relative to other phases or compounds. Since (/>'" of the listed ox- 
ides does not shift with pH value, when pH=7, the (j)"-" of ZnO and 
FeTi03 fall below (^(02/H20), indicating that these compounds are 
stable under illumination in the neutral or alkaline solution, which 
is consistent with the observation of ZnO in solution: no photocor- 
rosion at pH=10 while complete decomposition at acid pH=4.5— 

(ii) (p"-'' of all non-oxide semiconductors are higher than 



0(O2/H2O), at least 0.5 V more at pH=0, and the alignment of 
(j)"'^ relative to (^(02/H20) is not changed by pH value. This in- 
dicates these non-oxide semiconductors are thermodynamically 
unstable in aqueous solution and will be oxidized by the holes 
under illumination. The reason is also simple: the anions such as 
N^^, P-'^, As^^, S^^, Se^^, Te^^ can all be oxidized to neutral 
or positive valence states, e.g. N^^ to N2 or NO3", and S"^ to S 
or SO4 . It is interesting to note that the oxynitride TaON, which 
is between the oxide and nitride and naively expected to combine 
the better stability of the oxides and higher valence band (hence 
smaller band gap) of the nitrides, '^•29,30 j^^g jj^g oxidation potential 
close to the nitride Ta2N5, therefore it does not inherit the good ox- 
idative stability of the oxides. This further indicates that although 
doping or alloying the oxide semiconductors with weaker electron- 
negative anions can decrease the band gap, they become less stable 
with respect to oxidation in water. 

(iii) For the reduction potential (j)'''', all the non-oxide semicon- 
ductors have higher 0''' than (^(H+/H2) at pH=0, and thus are re- 
sistant to the electron reduction under illumination, corresponding 
to the situation in|2jb) with (j)'^" above (^(H+/H2). Since (^"° and 
(^(H+/H2) change at the same rate with pH value, this situation 
is not influenced by the pH change. Our results are consistent 
with the reported Pourbaix diagrams for CdS, CdSe, CdTe, GaP 
and GaAs,'4 which clearly show that their reduction and oxidation 
potentials are negative (higher in Fig.3) relative to ^(H+/H2) and 
(|)(02/H20) respectively at 0<pH<14. 

(iv) For (j)''" of the metal oxides, a few metal oxides have lower 
^'^ than (|)(H+/H2) at pH=0, such as CU2O, PbO, CUWO4, BiFeOs, 
BiV04, WO3 and C03O4, caused by the ease of reactions such as. 



Cu20 + H2^2Cu + H20 



(16) 



W03+H2^W02+H20 (17) 

In these compounds, the bonding between the metal cations and 
oxygen anions are weaker than the H-0 bonding in H2O, so the 
above reactions are energetically favorable. According to the ex- 
perimental Pourbaix diagrams,eli the WO3AVO2 and CU2O/CU re- 
duction potentials are V and 0.45 V relative to (^(H+/H2) re- 
spectively, in good agreement with our calculated (p'''^ . Although 
PbO, CUWO4, BiFe03 and WO3 have lower (j)'" than 0(H+/H2), 
they are still stable because their 0'"'' is above their CBM, and thus 
the conduction band electrons can reduce neither the water nor the 
semiconductors (Note (j)'''', CBM and (|)(H+/H2) change at the same 
rate with pH, so the stability analysis of these metal oxides is not 
changed by pH). However, CU2O, BiV04 and C03O4, have (j)'" 
lower than both 0(H+/H2) and CBM, corresponding to the align- 
ment inOb) with (/)'''' below (^(H+/H2), thus they will be reduced 
in the solution under illumination. Actually the reduction of photo- 
catalytic p-type C03O4 to inactive CoO had been observed exper- 
imentally.^' All the other oxides listed in [3] have (J)''" higher than 
(jyQi'^ fli2), and thus are stable against reduction. 

According to the above trends, one can find a series of metal ox- 
ides which are both resistant to the hole oxidation and the electron 
reduction and thus stable in the solution. Due to the lower CBM 
than 0(H+/H2), most oxides except CU2O, BiV04 and C03O4 can 
not be used as the photocathode, but can be used as the photoanode 
given their low VBM. Furthermore, the low VBM also makes the p- 
type doping difficult and most metal oxides are intrinsically n-type 
according to the doping-limit rule,— which states that a semicon- 
ductor is difficult to be doped p-type if its valence band is too low 
and is difficult to be doped n-type if its conduction band is too high. 
Therefore the metal oxides with low valence band should be better 
used as the n-type photoanodes. 

For the p-type photocathode, one can easily find metal non- 
oxides which have higher CBM than 0(H+/H2) and are also re- 



sistant to the electron reduction, thus can be used the photocathode 
if doped to p-type. According to the doping limit rule, the p-type 
doping in the non-oxides are easier than in the oxides due to the 
higher VBM. One may worry about the bad stability of the non- 
oxides with respect to the hole oxidation. However, we should note 
that, for the p-type photocathode, (i) the downward band bending at 
the semiconductor/water interface prevent the majority holes from 
reacting with water; (ii) the photo generated hole is expected to flow 
to the connected anode quickly through majority carrier conduction 
(at least when the device is working), - thus in this sense we can 
ignore whether it is oxidization resistant. As a result, some of the 
III-V and II-V semiconductors with suitable band gaps can be used 
as photocathodes, despite the fact they might be prone to oxidiza- 
tion from the pure thermodynamic point of view. The same can be 
said for n-type photoanodes where we are mostly concerned about 
their resistance to oxidization. 

Finally, we want to mention that the results represented in Fig.3 
considered only the thermodynamic resistance to the reductive and 
oxidative decomposition, but the real decomposition also depends 
on the specific kinetic processes. If the material is stable thermody- 
namically against the decomposition process, this material is stable 
regardless of the kinetic process. However, if the material is unsta- 
ble thermodynamically, there might still be ways to make it stable 
kinetically. For example, an oxide layer, which forms as a result of 
oxidation, but stops itself after a certain thickness, can be a good 
protective layer and make the material stable kinetically. For some 
p-type photocathodes, such protective oxide layer can also serve 
as a hole blocker, which prevents the recombination of the minor- 
ity electron and the majority hole, and thus can be essential to the 
efficiency of a photochemical cell. 

Note, through out our study, we have not considered the possibil- 
ity of the semiconductor compound to be dissolved in the solvent 
in dark condition. That can correspond to a reaction where the va- 
lences of cations and anions in the compound are not changed, but 
their total Gibbs energy is smaller in the solvation form than as in 
a solid crystal structure. For example, this can happen to GaAs 
within some pH range.— However, we do find some interesting 
cases for some compounds where reactions in the dark will occur. 
Comparing the oxidation potential (j)"-'' and reduction potential (j)''" 
in [3] most compounds have (j)"'' lower than (j)''", which is easy to 
understand since the oxidation means the electrons are taken away 
from the low valence band while the reduction means the electrons 
are added to the high conduction band. However, several excep- 
tions exist, i.e. CU2O, AlP, GaP, AlAs, Si and SiC have (j)'''' lower 
than (j)"-''. As a result, the electron-hole pair needed for the reduction 
and oxidization reactions can be generated spontaneously (they cost 
negative energy), without the help of any photon. This will cause 
these materials to decompose even in the dark. For example, the 
following two reactions are used to calculate the 0''' and (j)'"' of 
AlP respectively. 



2AtP + 3H2 -^2AI + 2PH3 (reduction) 



(18) 



2A//' + 3//20-s-A/203+2/' + 3//2 (oxidization) (19) 

giving that (p'''^ is lower than 0™, which makes the two reactions 
spontaneous, corresponding to a sum reaction. 



4AIP + 3H2O -^ 2AI+AI2O3+2P + 2PH2 



(20) 



The calculated Gibbs free energy change of this reaction is negative, 
indicating the reaction is exothermal, and thus AlP will be decom- 
posed in water even without any illumination. But once again, a 
kinetic barrier (e.g., an oxide layer) can block the otherwise ther- 
modynamically plausible reaction. 

In conclusion, we studied the thermodynamic oxidation and re- 
duction potentials for a series of photocatalytic semiconductors. 



and plotted their alignment relative to the valence and conduction 
band edges as well as the water redox potentials. According to these 
potentials, we found only metal oxides can be thermodynamically 
stable when used as the n-type photoanodes, and all the non-oxides 
are unstable due to easy oxidation by photo generated holes. This 
indicates although the oxides doped or alloyed by less electroneg- 
ative anions can have smaller band gaps, their stability is also de- 
graded. On the other hand, many non-oxides are resistant to the 
electron reduction, and thus may be used as p-type photocathode 
provided the hole oxidation is prevented in the working devices. 
The method we used is universal and can be applied to evaluate the 
stability of any semiconductors. We believe the presented stability 
diagram (|3} will be useful in guiding the search for stable photo- 
cathode and photoanode materials. 
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